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Abstract 

In this paper we present an analysis of the multiple hadron pro- 
duction induced by primary heavy quarks in e^e~ annihilation with 
the account of most complete and corrected experimental data. In 
the framework of perturbative QCD, new theoretical bounds on the 
asymptotically constant differences of the multiplicities in processes 
with light and heavy quarks are given. 

1 Introduction 

As it is already known in classical theory, the more is the mass of a charged 
particle, the less intensive is the radiation from it. In quantum field theory, 
in particular, QCD, this circumstance leads to a number of impressive ef- 
fects caused by heavy quarks, as, for instance, the effect of leading mesons 
which contain c- and 6-quarks in e^e~ annihilation [Ij. If we assume that 
the "radiation" in QCD is the radiation of gluons, which "materialize" af- 
terwards in form of detected hadrons, it is natural to expect, that in events 
induced by heavy quarks, the hadron multiplicity (except for decay products 
of the leading quarks) is smaller that in analogous events triggered by light 
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quarks. So natural may be an expectation that the difference between these 
multiphcities has to disappear at high enough energies. 

However, it is not always the case. For instance, even in QED, the dif- 
ference between cross sections (times the cm. energy squared) of fermion- 
antifermion pairs of different masses produced by two photons, does not 
disappear with increase of energy, but tends to a constant which depends 
on the masses of the final fermions [2] . The measurements of hadron multi- 
plicities in e^e~ annihilation associated with the primary quarks of definite 
flavors (in practice, u, d and s quarks are assumed to be massless) were car- 
ried out at lepton colliders with the collision energy of 29 GeV and more, in 
particular, at the SLC, KEK, LEPI and LEPII. 

In view of a great interest caused by these experiments (see below), and 
due to the presence of competing theoretical predictions, there is no doubts 
that similar measurements will be also done at future linear colliders, such 
as CLIC or ILC [3J. The first attempt to look for such effects was done in 
the framework of so-called "naive model". The essence of this model lies in 
the fact that an universal mechanism of multiple production of hadrons from 
some gluon system is adopted, which is insensitive to the quark flavors, while 
all the difference between processes induced by distinct quark-antiquark pairs 
arises due to a difference in energy available for the hadron production. 

Later on, another approach was accepted appealing directly to the cal- 
culational kit of perturbative QCD. However, the estimates obtained this 
time overestimated the data significantly. For all that, there was a clear in- 
dication of the asymptotic constancy of the multiplicity differences from the 
heavy and light quarks. Not long after, the same quantity, the multiplicity 
difference in the events with the light and heavy quarks, was calculated more 
accurately. The results appeared to be strikingly close to the experimental 
data. 

As was already said, the "naive model" [U [5] was the first attempt to take 
into account hadrons produced in addition to decay products of the heavy 
quark-antiquark pair in e~^e~ annihilation. Later on, it was argued [6] that 
the difference between multiplicities in events induced by heavy {Q = c, b) 
and light (/ = u, d, s) quarks. 
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Here (and below) it is assumed that we deal with mean multiplicities of 
charged hadrons, and "e" is the base of the natural logarithm (Ine = 1). 

The comparison with the data has shown that the "naive model" describes 
the data on 5bi up to Q = 58 GeV IH E El [9] , but underestimates the LEP 
and SLAG data ITTj [T2] . As for the formula which was obtained at the 
basis of the co-called "modified logarithmic approximation (MLLA)" ([2]), it 
has significantly overestimated both low and high-energy data on 5u- 

Detailed QCD calculations of the difference between associated multiplic- 
ities of charged hadron in e"'"e~ annihilation were made in [13]. The QCD 
expressions for 5qi from Ref. [13] appeared to be in a good agreement with 
experimental measurements of associated hadron multiplicities in e^e~ an- 
nihilation (see, for instance, [I11II5]). 

As we will see below, it is the hadron multiplicity in the light quark events 
that enables one to calculate the multiplicity differences 5qi. Recently, the 
data on average charged multiplicities in II events at different energies cor- 
rected for detector effects as well as for initial state radiation were presented 
in [THj. The corrected multiplicity differences averaged over all presently 
published results can be found in Ref. [15] : 

5^7^ = 3.12 ±0.14 , (3) 
= 1.0 ± 0.4 . (4) 

Because of the appearance of these corrected experimental data, a natural 
necessity has arisen to reconsider our predictions for 5u with the account of 
the data on the hadron multiplicity in the light quark events Nii{Q'^) as well. 

In Section [2] we define all the relevant quantities, since rigorous definitions 
are not usually presented by other authors. The analytical formula for the 
hadron multiplicity in e~^e~ annihilation is obtained. The QCD expression 
for the multiplicity difference is derived in the next section. In Section HI 
the upper and lower bounds on 6bi are calculated. In Appendix A the de- 
tailed derivation of the evolution equation for the multiplicity in the gluon 
jet is given. In Appendix B we discuss the connection between our approach 
and the scheme which uses the concept of the Altarelli-Parisi decay func- 
tions. The problems related to the gauge invariance guarantee which appear 
in perturbative calculations of the "light multiplicity" are considered in Ap- 
pendix C. 
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2 Hadron multiplicities in e^e annihilation 

The average multiplicity of hadrons in a qq event in the process of e^e 
annihilation is of the form 



where d'^^{k) is the propagator of the gluon with momentum k. Here and 
below (a, h) and (a', h') denote color indices. 

The subscript q denotes the type of primary quarks. In what follows, 
the notation q = Q {heavy quark) will mean charm or beauty quark, while 
the notation q = I {light quark) designates a pair of u, d or s-quarks which 
are assumed to be massless. In particular, Nii{Q'^) means the multiplicity of 
hadrons in light quark events, while Nqq{Q'^) is the multiplicity of hadrons 
in events when the process is induced by the heavy quark and antiquark of 
the type Q. 

The first term in the r.h.s. of Eq. ([5]), 2nq, is the multiplicity from the 
fragmentation of the leading quark (antiquark). It is taken from the analysis 
of data (2nc = 5.2, 2nb = 11.0 [B], and 2n; = 2.4 [H]). The quantity 
n^^ {k) in the integrand of Eq. is given by the diagram in Fig. [1] in which 
both the integration in the momentum of the final hadron and averaging 
in its polarization are assumed. This diagram corresponds to the following 
analytical expression: 



Here Jh{x) is the source operator of hadron h (spin indices are omitted), and 
/"(x) is the gluon (color) current [T7], 







x(0|(r/:'(x)J+(y))(T/J'(0)A(z))|0) 
= -^ [[[ d^xd'^yd^ze'^''Dl{y - z) 





6S 



S* , 



(7) 



= I 



6A^{x) 
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Figure 1: The diagram which describes the average multiphcity of hadrons 
with the 4-momentum h (sohd hnes) in the gluon jet with the virtuahty k"^ 
(winding hnes). 

where A'^{x) denotes the gluon field. D^{x) is the Pauli- Jordan function 

m- 

D^(x) =z(0|/i(x)/i(0)|0) , (8) 

where h{x) is an asymptotic hadronic field operator^ 

Since /c°/^(a;) = 0, and the final hadrons are colorless particles, we get: 

= {-9o.p + k^kp) 6^'''ng{e) , (9) 

where dimensionless quantity ng{k'^) describes the average multiplicity of 
hadrons in the gluon jet with the virtuality k"^ . It is, of course, gauge invari- 
ant, and depends only on the virtuality k"^. 

It is useful to introduce the average multiphcity from the gluon jet whose 
virtuality varies up to k"^: 

N,{k\Ql)= / 4n,(p2). (10) 

Very often Ng{k'^) is erroneously called the average multiplicity of the gluon 
jet with fixed virtuality k"^. This meaning should be addressed to rigik"^) only. 

The infrared cut-off Qq which separates perturbative and non-perturbative 
regions was introduced in flTUl) . We will use the "conventional standard" value 
Qq = 1 GeV for our numerical estimates (see Section H]). 

^For half-spin hadrons, DJ^{x) is replaced by S^{x) = {id + nihiD^lx). 
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The quantity Ug^k"^) cannot be calculated perturbatively. It is usually 
assumed that the average hadron multiplicity is proportional to UgihP'^Ql), 
i.e. the average multiplicity of (off-shell) partons with the "mass" Qq (the 
so-called local parton-hadron duality): 

ng{e)=ng{k\Ql)K{Ql), (11) 

where K{Q^) is a phenomenological energy-independent factor. The QCD 
evolution equations for both ng{k'^,Ql) and Ng{k^,Ql) are derived in Ap- 
pendix A. Let us stress, however, that the main results of the present paper 
(see Sections [3], S]) do not depend on explicit form of the function ng{k'^). 
In Eq. ^ the first factor of the integrand is given by 

Il%{q,k) = {-gnK;,M^^)^ (12) 

where n^^.^^(g, /c) is the two-gluon irreducible part of the relevant disconti- 
nuity of the four-current correlation functioij^ 

X (o| (TJ7(x)/;(y)) (rjr(o)/^(^)) |o) . (13) 

Here Jp™(x) is the electromagnetic current. The color current/^ (x) was de- 
fined above ([7]). 

In the first order in the strong coupling constant, n^^(g, fc) is given by 
two diagrams presented in Fig. [2] and Fig. [3] normalized to the total e"^e~ 
rate. This quantity is gauge invariant: 

Fn^t(9,^) = o. (14) 

Let us define 

^^(«^ - (^(- *=|.) / ^mvrn^' 

x(-g'--')S^tnfM\k\gk), (16) 

where Cp = {N^ — l)/2Nc, and A^^c = 3 is the number of colors. Then the 
average multiplicity in e~^e~ annihilation looks like 

K,{Q\ Ql) = 2n, + C,r ^ E{Q\ e) NgiP) . (16) 



^Note that 11°^.^^, is proportional to S^'' 
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Figure 2: The inclusive distribution of the massive gluon jet (winding hne) 
with the virtuahty k"^. The wavy hne is the virtual photon with 4-momentum 
q. The solid lines are light quarks. The cut lines correspond to on-shell 
quarks. 

The term E{Q^, k'^) is the inclusive spectrum of the gluon jet with the vir- 
tuahty up to k"^ emitted by primary quarksif] 

For the light quark case, the explicit form of E^Q"^, k"^) was obtained in 
Ref . ITS] . In terms of the variable 



fc2 



(17) 



we get: 



1 3 + 7(T /. 1 ^ ^ 



E{a) = (1 + 2(T + 2(T^) In ' {1 - a) - a{l + a) (\n 

a 2 \ (7 



+ 4(7(1 + a)/(a) , (18) 



where 



/(a) = --Li2(l + a), (19) 



and Li2(-2) is the Euler dilogarithm. 

■^It was explained in detail in Ref. pTSj that one should not consider this mechanism of 
hadron production via gluon jets as due to "a single cascading gluon" , as some authors 
believe [16]. That E is an inclusive spectrum of the gluon jets is seen, e.g., from the fact 
that the average number of gluon jets / {dk'^/k'^) E{Q'^, fc^) > 1. 
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Figure 3: The interference diagram which also contributes to the inclusive 
distribution of gluon jets. 



Let us introduce new variables 

= In 

k 



77 = In (20) 



and 

as well as the notation 



>^ = ln^, (21) 



N,ie) = CF^^N,ie). (22) 
Then Eq. ( IT6|) can be rewritten as 

Y 

N,g{Y) =2n, + Jd7]Ng{Y-r]) E{r]) = 2n, + N,{Y) . (23) 


The physical meaning of the function 

Y 

Ng{Y) = J dr]Ng{Y - f]) E{r]) (24) 



in Eq. ( l23l) is the following. It describes the average number of hadrons 
produced from virtual gluon jets emitted by primary quark (antiquark) of 
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the type q. In other words, it is the multiphcity in qq event except for 
the multiphcity of the decay products of these quarks at the final stage of 
hadronization (the first term in fl23|) ). 

The function E{ri) is presented in Fig. |H It has the asymptotics 

E{r^%^^ = E-^y^{r^) = r^-^-, (25) 

The derivative of E{ri) is positive for all r], as one can see in Fig. HI It 




Figure 4: The function E[ri). 



3 Multiplicity difference in QCD 

Now let us consider the difference between multiplicities in events induced by 
the light and heavy quarks, 5q;, which is defined by Eq. ([1]). The following 
representation was found in Ref. [T3] : 

5gf° = 2(nQ-nO-AiVQ(Fj. (26) 

Here the new notation is introduced, 

A7VQ(r„) =N^-Nq= [ dyN^iY^ - y) AEgiy) , (27) 
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as well as variables 

7/ = In 

k 



2 



2/ = lnT?, (28) 



and 



2 

mf 



Ym = \n-^. (29) 

Vo 

An important result which was obtained in Ref. [13] is that the function 

AEq = E-Eq (30) 

depends only on the variable 

p = exp{-y) , (31) 

but not on the energy Q. The explicit form of AEq is the following: 

7 1 20 

AEgip) = (l-3p+-p2)ln- + p(7p-20)J(p) + — -[1- J(p)] 

z p P ^ 

+ 7p+^, (32) 

where 



^In -^), p>4 



J(p) = <[ 1 , p = 4 , (33) 

/^arctan(^), p<4. 



Since AEq{ii) has the asymptotics 



A^Q(y) - y exp(-|y|) , (34) 



the integral in Eq. (|27|) converges rapidly at y ^ — oo. The function AEgd/) 
is shown in Fig. [5l We find that 

AEQ{y)\y^^ = AE-y-{y)=y-^-. (35) 

Another important relation comes from (l25l) and (l35!l : 

AEq(|/ - 1) - ^(y) ^ ^ ln2 exp(-|//2) . (36) 

oo Z 
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AEq (y) 



In other words, 



Figure 5: The function AEQ{y). 



/\EQ{y)c^E{y + l) (37) 



at large y. 

If one puts AEq^y) = E{y + 1), then (neglecting the contribution from 
the region y < —1): 

ANQ = Nu{mle)-2ni. (38) 
Correspondingly, the approximate expression for 6qi is of the form: 

5^^r = 2r^Q - Naimle) = 5^,^^^ , (39) 

where Sq^I"^^ is the MLLA prediction for the multiplicity difference [6]. We 
would remind that the function Nii{Q) describes the hadron multiplicity in 
light quark events at the collision energy Q. 

However, expression (!37j) is very far from the exact one in the region 
y < as it is clearly seen in Fig. [6l That is why, there is a large difference 
between 5^lla qqj^ expression (ES]). 



*For the beauty case, one has < 3.2. 
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AEQ(Y-1)-E(y) 
1 

0.8 
. 5 
0.4 



0.2 




Figure 6: The difference AEQ[y — 1) — E{y) as a function of tlie variable y. 

4 Multiplicity difference: upper and lower 
bounds 

In tliis section we will show that the difference between Sq}"^^ (l39l) and 5gf ° 
is indeed numerically large, and also obtain both upper and lower bounds on 

hi- 
lt is convenient to represent expression for ANq (1271) in the form: 





-1 



dyNg{Ym + l-y)E{y) 
dyNg{Y^-y) AEgiy) 



— oo 

Ym + 1 



+ / dyNg{Y^ + l-y)[AEQ{y-l)-E{y)] 



= [Nairn' e) -2ni] + 5iV« + 6N''> (Y^) 



r(2)/ 



(40) 
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that results in the formula (see Eq. 

= 2nQ-Naimle)-6Ni^\Yj-5N^^\Y^) 

= Si;r^-K\Yrn)-SN§\Yj. (41) 

Here we have introduced notations 

-1 

N^Qiym)= jdyN,iY^-y)AEQiy), (42) 



and 



N;^\yJ= j dyNg{Y„.+ l-y)[AEQ{y-l)-E{y)]. (43) 







Note that both NQ^(Ym) and n'q' (Ym) are positive since AEQ{y) > at all 
y and AEgi^y - 1) - E{y) > at ^ (see Fig. and Fig. E>- 
In order to exploit the corrected data on Nif^Q"^) at Q = 8 GeV, 

NiiiS.O GeV) = 6.70 ± 0.34 , (44) 

we assume that = 4.85 GeV, that corresponds to mbi/e = 8 GeV. 

The estimates show that the dominant correction to Sgf^ is SNq \ not 

SNq\ To calculate a lower bound for SNj^^\ we will use the following in- 
equality in the region y ^ 0: 

AEQ{y)^E{y + AyQ). (45) 

The quantity Ayq is a solution of the equation 

AEQ{X^)=E{Y^ + AyQ), (46) 

where Y^ is defined above (1291) . Then we get from Eqs. ( H3l) and ( l45l) : 

Aj/Q-l 

6N^^^-^Na{Y^ + AyQ)~Na{Ym + l)- j dy Ng{Y^ + AyQ-y)E{y) . {Al) 



For our further estimates, we need to know the hadron multiplicity in the 
light quark events in the energy interval 2.5 GeV ^ Q ^ 28 GeV. By fitting 
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the data on hadron multiplicity in the light quark events at low Q, we get 
the expression: 

Ar,;-(Q2) = 2.07+ 1.11 In Q + 0.54 In^ Q . (48) 

Putting Qo = 1 GeV, we find Ayi, = 1.61. Taking into account that the 
last term in Eq. ( l47ll is negligible]^ we get from (l47l) . (l48l) : 

5nI^^ ^ 1.07 . (49) 

Correspondingly, our prediction accounting the revised data on the multi- 
plicity in events induced by the light quarks, 

5,^^° ^ 2nb - NaiYm + Ay^) = 3.33 ± 0.38 , (50) 

appears to be lower than our previous result 6m = 3.68 [13]. We used the 
experimental value 

2rab = 11.10 ±0.18 . (51) 

The error in measurements of Nij was taken to be ±0.34. Let us stress 
that our upper bound (l50l) is very close to the present experimental value of 

Sir m- 

Now let us derive a lower bound on 5^*"°. To do this, let us start from 
Eq. (127|) . It is convenient to represent the integral in (1271) as a sum of two 
terms 

AN, = JdyNg{Y,-y)AE,{y) + JdyNg{Yrn-y)AE,{y) 

-4 -1 

= AiV(^) + AiVf , (52) 

where Yf, = ln(m^/(5o) — 3.16. Consider the first term in (1321) . One can 
check that 

AE{y) < 0.18 E{y + 5.8) (53) 
in the region —4 < y < —1, that leads to the inequality 

4.8 

AN^'^ < 0.18 J dyNg{Y, + 5.8 - y) AE,{y) . (54) 

1.8 



'^Since E{y) < 0.02 in the region ^ y ^ Ay^ - 1 = 0.61. 

®We take into account that the region — oo < y < —4 gives a negUgible contribution to 



ANh 
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The estimations show that Ng{Yiy + 5.8 — y) < 2 iVg(4.8 — y) when y varies 
from 1.8 to 4.8. Thus, we get: 

AN^'^^ < 0.36 [NiiiQ = 11 GeV) - NaiQ = 2.5 GeV)] = 1.54 ± 0.17 . (55) 



The second term in (1521) can be estimated by using the inequahty 

AE{y) < 0.62 E{y + 3.5) (56) 
which is vahd in the region —l<y<Yb. Then 
AA/'f ^ < 0.62 [NiiiQ = 28 GeV) - NaiQ = 3.5 GeV)] = 4.61 ± 0.30 . (57) 



As a result, we obtain from Eqs. (1261) . (!27j) and (1551) . (1571) the lower bound 



on 6, 



QCD 
bl 



> 2.55 ± 0.39 



(58) 



Fig. [7| demonstrates that our QCD predictions are very close to the average 



measurement Slf^ = 3.12. 
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Figure 7: Our QCD result for 6hi (a corridor between the two solid lines), 
and the MLLA prediction [16] (a corridor between the two dashed lines) vs. 
experimental data. 
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Our results can be compared with the recently published new MLLA 
results [TB] : 

The next-to-MLLA results are: 

^NMLLA ^ 2.6 ± 0.4 , (60) 

Thus, the lowest-order MLLA expression ([2]) is not stable against higher 
order corrections. As it was said above, the formula which was used in the 
MLLA approach fl5^ is only a part of the exact QCD formula ([T]), (1271) 
in the approximation E{y) = AEgly — 1). This approximation is quite 
a rough one (see Fig. [H]), and the deviation of the function AEqly) from 
AEQ^"^{y) = y-3/2, as well as the deviation ofE{y) from E^^y'^iy) = y-1/2, 
cannot be neglected. 

Note that the limit y oo means that the invariant mass of the gluon 
jet k"^ tends to zero, since k'^ = Q^exp(— (see f l25]) ). We also think that 
using the argument k\ in A^^^ instead of k"^ in the MLLA scheme (see [TB] 
and references therein) has no reasons. Indeed, is a time-like vector (/c^ = 
/cq — /c^ — > 0), and the gluon jet with large transverse momentum k± has 
the small invariant mass k'^. In such a jet, the multiplicity is small, since the 
phase space for secondary particles is actually defined by the invariant mass. 



5 Conclusions 

The formula for the difference between hadron multiplicities in e~^e~ annihi- 
lation into hadrons induced by light and heavy primary quarks (with Q is a 
type of a heavy quark) is derived: 

= 2nQ-Nuimle) (62) 

-/;^f...)[A.«(f)-.(t)] 

Here Ng{k'^) describes the average number of charged hadrons in the gluon jet 
with the virtuality which varies up to k"^, and E, AEq are known functions. 
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By using the data on the hadron multiphcity in hght quark events Nii, 
corrected for the detector effects and initial state radiation effects [15], we 
have obtained from fl62l) the bounds: 



2.2 < 5^^^° < 3.7 . (63) 

We would like to emphasize that our estimate does not depend on a specific 
choice of the function Ngik"^), and it is in a good agreement with the average 
experimental value = 3.12±0.14. 

Two last terms in (1^21) . which are subtracted from the first one, are 
positive and numerically large. In particular, for the case of the beauty 
{rriQ = TUb, uq = Uh) the second term in (l62l) (dominating the third one) is 
equal to 1.1. 

As a result, the deviation of the MLLA prediction, 

5,^^^^ = 2n, - Naimle) , (64) 
from the exact expression^ 

Sr = 2(n. - n,) - I^^^N.ik^) ^E, , (65) 

appears to be quite significant. 

Let us also mention that our numerical predictions for the charm quark 
case, 

5j^^{Q = 170 GeV) = 0.99 , (66) 



were derived in [T3j before the precise measurements of 6ci were made [TT] 
As one can see, our value (166|) is in a nice agreement with the average exper- 
imental value (jlj) (see also Fig. [7]). 

Some part of the results of this paper was published in [18]. We are 
indebted to W. Ochs, the correspondence with whom stimulated, to some 
extent, the appearance of the present paper. 

The work is partially supported by the grant RFBR-06-02-16031. 



^This formula is an equivalent compact form of Eq. for Q = b. 

the low energy measurements HllBl, the total error of 6ci was about ±1.5. 
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Appendix A 



Let us consider the average multiplicity of off-shell quanta with the "mass" 
Qq in a gluon jet whose invariant mass is p^. It obeys the following integral 
equation [TOll^: 

JQI I- Jo 

X ^^P,,{z)n,{l\Ql). (A.l) 

Here k"^ is the virtuality of the parent quark which emits this gluon jet, and 
Pgg{z) is Altarelli-Parisi time-like decay function. 
The inequality 

in Eq. (lA.ip is a kinematical bound, while the bound 

z ^ 1^ (A.3) 

is a dynamical one. The latter is nothing but the angle ordering rewritten in 
terms of momentum fractions and virtualities (see Refs. m, M)E 

Assuming that parton shower develops mainly via soft gluons, one can 
put in flAAjl 

2Ca - . (A.4) 



Then the sequence of equations [21] , 
N,{k\Ql)-l 



dp'^ f^^ dp r dz 2 i2\n/ i2 2\'^s{zp'^^ 



I ^ I ^ I -e{zp'-ty{ze-p')'^^^nS\Ql) 

JQI P JQl ^ Jo z TT 
jQl P ^ JQl I Jo Z TT 



^Namely, the emission angle of the secondary gluon with the virtuality P is less than 
the angle of the primary gluon emission off the parent quark. 
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JQI P ^ [JqI I Jo ^ 

+ / ^e{kH'-p') / -e{zp'-i 

Jo?, I Jo ^ 



Ois{zp^) 



^gi} 1 Qo 



,2 TT 



I P 



dr"^ 







df 



2i2\ 



'^eir'e-p^) I ^ :^e{p^-kH 



dr'^ 



2 



df 



^eir'k'-p") / ^eikn'-p"") 



2 p 





TV 



JQl P JQ 



^ 



dp" 



+ Ca -^e{p'-Jk^Ql) 9ik'r'-p 



2 

^ 

dp" P'dr^asir^] 



dr'^ 



2 





TT 



2 



2 TT 





Q2 7q2 r2 27r 



results in the following formula for Ng{k , Qq): 



Ng{k\Qi) = l + CA j % 

^ 



.2 2tt 



From (1A.6I1 we obtain the differential equation: 

(^'^f) ' = CA^Ng{k\ Ql) 



with the boundary conditions 

Ng{k\Ql 



k'^ — Qq 



1, k'^^Ng{k\Ql) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 
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This equation has the solution: 



N,{k\Qt) 



2Ca, k- 



2Ca , Qo 



+ h 



— — In — \ Kn 
nb A? 



with the asymptotics 



2Ca , Ql 



2CU ^2 
nb ''A2 



(A.9) 



(A.IO) 



Here b — (33 — 2A^j)/127r, where Nf is a number of flavors. 

Let us stress that the equation for an isolated gluon jet (which has no 
parent parton with virtuahty k"^) would be of the form: 



X 



that leads to the formula: 





'-P,,{z)nS\Ql), 



27r 



(A.ll) 



N,{k\Ql) = l + 2CA r^^^N.ir^Ql). (A.12) 



This equation results in a wrong expression which does not take into account 
the interference effects: 

^.(^^Qo)L.»Qg - exp {2^^^ . (A.13) 



Appendix B 

Now we will reproduce the asymptotic relation between the average multi- 
plicity in the light quark event and that in the gluon jet by using the Altarelli- 
Parisi decay functions. Let / be a 4-momentum of the primary quark which 
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emits a massive gluon jet. The ladder diagram in Fig. [2]leads to the equation 



X 



2 



^.(fc^Qo) (B.i) 



(for simphcity, here and below we omit the contribution from the leading 



hadrons, 2nq) 



In the leading logarithm approximation, 

one comes to the expression (r^ = z/^): 

JQI ^ JQl ^ 

)2 vr L2 r 







2 



'^^ / ^:2 In — A/^5(r ,Qo) , (B.3) 



2 TT 





where we have used the relation: 



NAr^Ql)= f%n,{k\Ql). (B.4) 
•JQl P 

The integral equation for Ng{r^,Ql) has been obtained in Appendix A (see 
Eq. dAJD). 

The formula (1B.3P can be represented a^ 

NniQ', Ql) =CFr ^ ^ E-y-{Q^ e) N,{k\ Ql) , (B.5) 

with the function E^^^'^{Q'^ , k"^) defined above fl25|) . It is worth to compare 
this approximate expression with our exact formula f|T6l) . 



We have added the non-leading term -1/2 to ln{Q /r ) in deriving (jB.sp from (jB.Sp . 
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At large Q^, we obtain from (1B.3P and (lA.lOp the well-known asymptotic 
relation: 

Nu{Q\Ql) ^ ,^^N,{Q\Ql). (B.6) 

Remember that Ng{Q'^, Q^) describes the average number of virtual partons 
in the gluon jet whose invariant mass varies from Qq up to Q. 



Our formula (11611 should be reproduced from the starting equation (IB.ip 
provided: 

1. Contributions from both the ladder (Fig. [2]) and interference diagrams 
(Fig. [3]) are added together in deriving the expression for Nij^Q'^, Ql) 

2. Both non-singular terms in variable z and power corrections 0(A;^//^) 
are taken into account in Pqg{z, k'^/P) 



Appendix C 



The average number of hadrons in e~^e~ annihilation is, of course, gauge in- 
variant quantity. However, in perturbative QCD we calculate the multiplicity 
of virtual partons. Apriori one can not be sure that it does not depend on a 
gauge. 

Unfortunately, this important problem was not studied to a considerable 
extent by other authors. As a few exceptions, let us mention Refs. [22] 
and [23] • In the first paper the gauge invariance of the multiplicity in e+e" 
annihilation has been checked in one-loop approximation. In the second 
one the gauge dependence was considered for a fixed coupling constant and 
without accounting for interference effects. 

Below we will analyze a possible gauge dependence of the partonic mul- 
tiplicity in the case of light primary quark in a class of axial gauges. It is 
given by Eq. fIB.ip in the gauge n^A'" = with the gauge vector 



= ^(1,0,-1) (C.l) 

(z-axis is chosen along a 3-momentum of a primary quark in the c.m.s. of 
colliding leptons). The argument of the decay function Pqg{z) in (IB. II) is the 
ratio 

. = ^^. (C.2) 

In 
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The emission of the gluon jets from the primary antiquark is suppressed in 
this gauge (IC.ip . 

Analogously, in the gauge 

n^ = -^(1,0,1) (C.3) 

the massive gluon jet is emitted by the antiquark while its emission from the 
quark is suppressed. 

Let us choose the gauge in which both quark and antiquark make com- 
parable contributions to the emission: 

= (1,0,0). (C.4) 

Then Eq. ( ]B.1|) is modified as follows: 

Nu{Q\Ql) = -w dz-^'^' 'I - u .... n.2 ^2. 



2 P 



^2 .,2 /•! 



2 r Jri2 



1 
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+ I ^ ^ I ^ dz^^ I ^P,,{z)n,{k\Ql), (C.5) 



where 

P„U) = P,,(^) = 2Ca — (C.6) 

at small z. These relations mean that the massive jets are emitted by the 
quark and antiquark with the same probability in the gauge (lC.4p . 

By omitting non-leading terms, we get from (IC.5I) . flC.6P the expression. 



2C, [ dr'^^N,{r\Ql) [ 



f^'dr'asir') Q\ ^ ^2^ ,p 7, 

^ 'q2 i^7^^5(^ ^Qo) , (C.7) 
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which coincides with formula (IB. 30 . 

Note that in the gauge (lC.4p the decay function in evolution equation (lA.ip 
is of the form: 

~ 2C^ . (C.8) 

This results in the effective cut on the integration variable z from below: 

> ^ . (C.9) 



As one can see from (lA.SP and (]C.9|) . it is the dynamical bound (1A.3|) that 



smoothes out the singularity of Pgg{z) at z = 0, but not the bound (]C.9P 
arising from the gauge vector. The latter can be safely omitted in Eq. (lA.ll) . 

Thus, we conclude that both the relation of the light quark multiplicity 
Nii{k'^,Ql) with the gluon multiphcity Ng{Q^,Ql), and the evolution equa- 
tion for Ng{k^,Ql) do not depend on the gauge vector n^. This conclusion is 
also valid for a general case, = (rio, 0, n,,), where uq ^ in,, [191 [20]. Note 
that the proof of the gauge invariance needs an account of the destructive 
interference in the emission of the gluon jets, that leads to the condition 
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